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I. INTRODUCTION 

The Frenkel-Kontorova model of classical particles connected to their nearest neighbors by springs and 
placed in a one-dimensional periodic potential V [Q-D is of interest for several different reasons. First, it 
provides a simple example which exhibits some of the complex phenomena (discommensurations, devil's 
staircase phase diagrams) associated with commensurate-incommensurate phase transitions For some 
recent studies, see [^,0; also ^ for extensions to quantum systems, and |^ for classical particles at a finite 
temperature. Next, various interesting dynamical phenomena result if forces, either static or dynamic, 
are applied to the particles in a Frenkel-Kontorova chain; for some recent work, see 

[|0|I3- Finally, the 

equilibrium configurations of this model form the orbits of an area-preserving twist map P4] , so that their 
study is part of the subject of nonlinear dynamical systems p3[ |. 

A special, exactly solvable realization introduced by Aubry Q and Percival |l^, in which V consists 
of a potential made up of identical parabolas, has played a pioneering role in the study of general Frenkel- 
Kontorova models, ft is particularly useful for understanding the properties of incommensurate ground states 
in the "Cantorus" regime: such configurations define invariant Cantor sets of the associated area-preserving 
map, on which the motion is quasiperiodic with an irrational winding number w. Cantori, or Aubry-Mather 
sets ll^, have been the subject of considerable research ||l7|- ^ , as they constitute the ubiquitous remnants 
of broken KAM tori in non-integrable Hamiltonian systems | |23[ . 

The present paper is dedicated to a generalization of this solvable model to the case in which the parabolas 
all have the same curvature, but need not be uniformly spaced nor located at the same height. The only 
constraint is that after N parabolas, where N is finite, the whole pattern is repeated, thus giving rise to a 
periodic V . See the example with iV = 3 sketched in Fig. |l|. The case iV = 1 is, of course, the one studied 
earlier by Aubry and Percival. 

The generalization is important for two reasons: First, as we pointed out for iV = 2 in a previous 
paper p4| , qualitatively novel phenomena emerge which are not present in the case iV = 1. Examples are 
incommensurate defects, discontinuous Cantorus-Cantorus phase transitions, and independent orbits of gaps 
composing the complement of the Aubry-Mather set. Second, it provides a tool, through the possibility of 
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approximating a smooth potential by a set of parabolas, for identifying true ground states in the presence of 
infinite sets of competing metastable Cantorus configurations, as encountered in Frenkel-Kontorova models 
with smooth V differing from a simple cosine . 

In the following sections we give the detailed solution for general N , which in essence consists in replacing 
the original infinite-variable minimization problem by that of finding the minimum of a certain convex 
function defined on a finite simplex. While the second minimization problem is not altogether trivial, it can 
be carried out to quite adequate precision in particular cases by a combination of analytical and numerical 
techniques. Efficient numerical methods are the subject of the paper which follows. 

Previous research on this topic of which we are aware includes an unpublished thesis of F. Vallet |^. He 
considered some particular cases for > 1, including the N = 2 situation with equally spaced parabolas 
discussed in our previous publication, for which he wrote down an explicit expression for the ground state 
energy. And he pointed out that the same methods could be applied to the general N situation using a 
parametrization of the hull function very similar to the one employed here. More recently Kao, Lee, and 
Tzeng [|7| have addressed the problem of finding the ground state by using a hull function to find solutions 
of the equilibrium equations. Their results are consistent with those in the present manuscript. 

In addition, Aubry, Axel and Vallet considered a somewhat different extension of the = 1 solvable 
model, with even and odd numbered particles subjected to different potentials While there are certain 
similarities, this model (and the method used to solve it) seems to be distinct from the N = 2 case considered 
here. Recently Kao, Lee, and Tzeng have solved a model in which the parabolas have negative curvature; 
the results are very different from the case of positive curvature discussed here. 

Our material is organized as follows. Section || contains basic definitions for the model including the 



explicit form of the potential V . In Sec. Ill the energy for an equilibrium state is rewritten with the help 
of a Green's function in a form which is combined with some of Aubry's exact results for ground states in 
Sec. 0. The result is used in Sec. ^ to set up an explicit minimization problem of a convex function over 
a finite simplex, the solution to which yields both the energy and configuration of the ground state. The 
conclusion. Sec. indicates a number of possible generalizations. Minor technical issues are discussed in 
two appendices. 

II. ENERGY 

We shall consider a system described by an energy 

^= E iliuj+i-u.f + Viu,)], (2.1) 

j=-oo 

where Uj is the position of the j-th particle, while 



V{u) = V{u + P) =mmp{n;u) (2.2) 

n 

is a periodic potential made up of segments of parabolas, as illustrated by Fig. |l|. The n-th parabola is given 
by the formula 

p{n;u) = ^k{u - tnf + hn- (2.3) 
We assume that the positions of the minima of these parabolas form a monotonic periodic sequence, 

tn < tn+1 ; i,i+JV = P + tn (2.4) 

and the values at the minima are periodic, 

hn+N = hn, (2.5) 

but otherwise arbitrary. 

Because all the parabolas have the same curvature, the u axis consists of segments separated by a sequence 
of dividing points, (see Fig. |l|), 

tn < tn+l] ^ri+AT = P + tn, (2-6) 
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with the property that for u in the open interval (which may be empty), the minimum in (2.2) 

is achieved for n = m, and for no other value of n. Thus this interval is uniquely associated with the m'th 
parabola, and if Uj fal ls i n this interval, we shall say that particle j is "in" parabola m, or n{j) — m. 
Given the periodicity (2.6), it is clear that two such intervals whose subscripts differ by N will be of the 
same length, and will have corresponding values of V{u). Thus we shall say that such an interval, or the 
corresponding parabola, is of "type l<l<N,iim = /(mod N). 
Let 

Tj = t,^(^j) (2.7) 

be the minimum of the parabola containing the j-th particle, in the sense just discussed. Then ( ^.l| ) can be 
rewritten as 

H = Ho + H,, (2.8) 

where 



j=-oo 

and 

oc 

^1= E '^"O)- (2.10) 

j = -oo 

For later analysis it is convenient to rewrite (|]^) in terms of 

^j=%-^i> (2-11) 
the distance of the j-th particle from the minimum of the parabola which contains it. One obtains 



oo oo oo ^ 



J — — oo k—~oo J — — oo J — — oo 



where 



and 



= ^j+i - Tj, AVj = Tj+i - 2tj + rj_i = Atj - Arj_i (2.13) 
i/a(0) = 2 + «. Hail) = Ha{-1) = -1, iJ^ (fc) = for | fc | > 1 . (2.14) 



III. EQUILIBRIUM STATES 



An equilibrium state is one in which the force is zero on every particle, which means dHQ/duj = for 
every j, and thus, by ( ^.12 ), 



J2 Ha{j - k)uk ^ 



(3.1) 



k— — c 



If we assume that the A^r, are known and are bounded as a function of j, the solution of (3.1) may be 



written as 



(3.2) 
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using the Green's function (decreasing at infinity) 

G(fc) = e-''l'=l/2sinhr, 

which satisfies the equation 



£ Ha{j -k)G{k-l)^5. 



k— — oc 

Here r is a positive number defined by: 



K = (2sinhr/2)^ 



Multiplying ( |3.fD by Uj and summing yields, see (2.f2) and (3.2): 

OC CX3 OO OO 



3 = -oo 



J — — OO /c — — OO 



J = -oo 

With the help of the function 

9{j) = G{J + 1) + G(j - 1) - 2G(j) + 5,0 = (tanhr/2)e-'-l^l 
one can rewrite (13.21) and (|3.6|) as: 



E 9U~k)Tk, 



k— — OQ 



OO OO 

^0=9 E E 9{3-k)ATjATk. 



j — — oo k— — oo 



The fact that 



E 9U) = 1 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



means that (3.8) expresses the position of particle j as a weighted average of the positions of the parabolas 
occupied by the different particles. 

It is useful for later purposes to extend g{j) to a function g(x) defined on the entire real axis by means of 
linear interpolation between successive integers, (see Fig. which is to say: 



]{x) := [1 — frac(a;)]5(Int(a;)) + frac(x)5(l + Int(a;)), 



(3.11) 



where Int(x) is the largest integer not greater than x, and frac(x) stands for x — Int(a;). Note that the 
resulting function is convex for a; > 0, and also for x < 0, and possesses a piecewise constant derivative with 
discontinuities whenever x is an integer. 

If the equilibrium configuration is monotonic in the sense that for all j , 



nU + 1) > n{j), 

that is, Atj > 0, then one can write Hq, (|3.9[) in an alternative form with the help of the function. 



maxjjeZ | n{j) < m}. 



(3.12) 



(3.13) 



the number of the last particle in parabola m, or of the last particle preceding this parabola if it is empty. 
(One can think of Im as a sort of inverse of the function n{j).) If 



(3.14) 
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is the distance between the minima of two successive parabolas, it is easy to show that 



(3.15) 



m— — oo 



Thus if particle j is not the last particle in some parabola, which is to say particle j + 1 is in the same 
parabola, /,„ never takes the value j and At, = 0, so both sides of this equation vanish. If j is the last 
particle in parabola m and parabola m + 1 is not empty, then Ar^ = At,„, = j, and Im+i is larger than 
J, so again (3.15) holds. If, on the other hand, parabola m is followed by one or more empty parabolas, then 
for each of these / takes the same value j, and thus Ar, is expressed as the correct sum of At intervals. 
With the help of ( p^ , (U) can be written as 



Ho 



E 



E 

m— — oo n— — oo 
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^ rn. 

)AtmAt 

r. 



(3.16) 



with sums over parabolas instead of particles. Similarly (2.10) may be written in the form 

oo 

Hl^ (Irn - Im-l)hm, (3.17) 

m— — oo 

because /,„ — /m-i is equal to the number of particles in parabola to. 



IV. CONTRIBUTION OF Hq TO THE GROUND STATE ENERGY 



For an infinite configuration {uj}j^z of particles, expressions such as (2.1 ) ar e formally infinite and do not 
possess a direct interpretation. Following Aubry |3y], we shall understand (2J.) as prescribing the change in 
energy if any finite set of particles are displaced from one set of positions to another. A minimum energy 
configuration is one for which any such change increases H (or leaves it constant) and a ground state is 
a minimum energy configuration which satisfies certain recurrence conditions pO| |. Aubry has shown that 
under suitable conditions (see also the additional remarks in App. ^ any minimum energy configuration 
has a well-defined average separation between particles. 



uj = Pw lim (uk — Uj)/{k — j), 

(fe-j)^oo 



(4.1) 



independent of how j and k behave individually (thus j = —k, or j = k/2 give the same result when k — > oo). 
We shall call w = lu/P the "winding number" , a notation clearly motivated by the twist map analogy jl^ . 
For a given uj, there is a well-defined energy per particle 



lim 



1 



K-l 



E[J( 



{K-J)~,oo (K - J) ^'2 



Uj+i ~Ujf + V{uj)], 



(4.2) 



where {uj} is any minimum energy configuration for this lo. (Note that this is diff erent from the process of 
choosing, for each J and K fixed, the {uj} which minimizes the finite sum in (4.2) ! ) 
If w is irrational, any ground state configuration with uj = Pw has the form 



f^{juj + c). 



(4.3) 



where c is some constant, and fui{z) is a monotone strictly increasing "hull" function ||30(| which is step 
periodic, i.e.. 



fJz + P)=P + U^. 



(4.4) 



When is discontinuous, as it is for the V(u) considered here, the discontinuities form a dense set, and 
there are two versions of f^^: /+, which at each discontinuity takes the maximum possible value consistent 
with monotonicity, and is thus right-continuous, or upper semi-continuous, and /j, which takes the minimal 



5 



possible value at each dis continuity, and is thus left continuous, or lower semi-continuous. Either /J" or /~ 
may be employed in (^^), though not both simultaneously [^ . 

When w = p/q is rational, with p a nd q relatively prime integers and q > 0, one can again write ground 
state configurations in the form ( [4.3| ), where the function foj{z) is defined on the discrete set of points 
z = Ps/q, where s is a ny in teger, and c is an integer multiple of P/q. Defined in this way, fu>{z) is strictly 
increasing and satisfies (4^). See App. |^ for the derivation of this result, which does not seem to be explicitly 
stated in For convenience of exposition, we shall assume that w is irrational in the following discussion, 
as the extension to the case of rational w is straightforward. 

In a ground state configuration no Uj can come closer than some fixed finite distance to one of the dividing 
points tn, ( |2.6| ) (see App. ^). Hence fuj{z) will have jump discontinuities with every t„ falling in the interior 
of one of these discontinuities, as illustrated in Fig. |3[ As fujiz) is monotone, this means that the z axis is 
cut up into segments by dividing points 



Zn ^ Zn+l\ Zn+N — Z„ + P 



with the property that 



Z < Zn < Z 



(4.5) 



(4.6) 



Thus there is a correspondence between intervals on the z axis and those on the u axis in that (fn_i, Zn) is 
mapped by into the interior of (t„_i,t„), so that intervals of type / on one axis correspond to intervals 
of type I on the other. 



The preceding discussion allows us to find an explicit formula for the integers {Im} defined in (3.13) in 
terms of the Zm- Assume the hull function is /~. Then, see Fig. |^, the number of the last particle falling in 
parabola m (or preceding this parabola if it is empty) is the largest integer j for which 



which means 



jcu + C < Zn 



I,n = Int[(z„ - c)/uj]. 



(4.7) 



(4.8) 



(If one employs /+ rather than , (4.7) is a strict inequality, and Int in (^) should be replaced by int, 
where int(a;) is the largest integer strictly less than x. The energy per particle is independent of the choice 
between /+ and .) 

The ground state energy per particle corresponding to Hq, ( [3.1(: ), can be computed as follows: Consider 
a set of LN successive parabolas, wher e L is large. Suppose that two integers toq, hq between 1 and N are 
given, and consider summands in (3.16) corresponding to pairs 



{m,n) = {mQ + fiN,no + {^i + iy)N), (4.9) 

where /i and v arc i nteg ers, and fi is between 1 and L. By ( ^.4| ), At„j and At„ are the same for all m and n 
of this type, while (4^), the step periodicity of the z„, allows us to write 



In -Im= l^t{D + A + ^/w) - Int(A + fi/w), 

where w is uj/P, 

A = Cma/w-c/uj, D ^ {Cno - Cmo + 

and we have introduced the notation 

Cn ■ Zn/P- 



(4.10) 



(4.11) 



(4.12) 



As ^ varies, the right side of (4.10) takes on only two values: 1 + Int(D), which occurs for a fraction 
frac(r') of the values of /z, assuming w is irrational and L is very large, and lnt{ D), which occurs for a 
fraction 1 — frac(-D) of the fj, values. Hence the sum of g(/„ — Im) over the pairs ( |4.9| ) with v fixed is for 
large L given by Lg{D), using the extended definition of g in ( [3.11 ). 
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What remains is a sum over v, which can be extended from — oo to oo with neghgible error because g cuts 
off exponentially, and sums over mg and ng. Upon dividing by the numbers of particles LP /to = L/w in a 
system of length LP, one obtains 



eo = ^ ^ ^ Ai„,A<„g(C„ - Cm) (4.13) 



■m— 1 n— 1 



for the energy per particle where the subscripts of (4.9) have been omitted from m and n, and we have 
introduced a new function 

oo 

g{x)^w g[^] (4.14) 

z^— — oo 

which depends implicitly on both w and k (or r). Note that 

g{x) ^g{i + x) = g{-x) (4.15) 

is periodic and symmetrical about x = and x = 0.5, and convex for x between and 1, as depicted in 
Fig. |. 



V. FINDING THE GROUND STATE 



The energy per particle ei corresponding to Hi, ( 2.10 ), is easily evaluated when w is irrational by n oting 
that the sequence of equally spaced points on the z axis which form the arguments of /^^ in (4.3) will 
eventually approach a uniform distribution if mapped modulo P into the interval [0, P). That is to say, the 
fraction of these points falling inside an interval of type I on the z axis, which is the same as the fraction of 
particles whose positions lie in parabolas of type Z, is equal to 



ipi = Ci- 0-1, 



with C/, ( [4.12D , equal to zi/P. In view of (|45|), 

Cn+N = 1 + Cn 



and therefore 



N 



1=1 



As a consequence. 



N 



N 



1=1 



1=1 



where 



Vi = hi+i - hi, 



and because the /I's are periodic, ( p.5| ) 



N 



(5.1) 
(5.2) 

(5.3) 



(5.4) 



(5.5) 



(5.6) 



1=1 



Thus the total energy per particle e = ep + ei can be written in the form 
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1 ^ 
-T 

2 ^ 

n=l 



N 



(At„)2g(0) 



N-l 

E 

m— 1 n— m+1 



AT 

E 

m— 1 



(5.7) 



by combining (4.13), slightly rewritten using (4.15), with (5.4). Here hi is an additive constant which is 
set equal to zero in the following discussion. Then e is a function of the unknown parameters CijC2j ■■■,Cn, 
and must be minimized by varying these parameters. In fact, adding the same constant to every C„ leaves e 
unchanged (note Eq. 5.6) and therefore, since the -0's, (5.1), cannot be negative, it suffices to consider the 
problem of minimizing e on the simplex 



- Co < Ci < C2 < ... < Cjv-1 <Cn = 1- 



(5.8) 



One can think of the {Cm} as the positions of a set of N "quasi-particles" located on a circle of unit 
circumference, interacting with each other through a set of pair potentials given by Q and, in addition, 
subjected to constant single-particle forces (the r/'s), whose sum is zero. The fact that G{x) has a minimum 
at a; = ^ means that the pair forces tend to keep the quasi-particles separated, although the imposition of 
suitable single-particle forces can push two or more together. However, the order ( ^.8[ ) must be preserved: 
these quasi-particles have hard cores and will not move through each other. 

), is a sum of continuous convex functions. 



Note that e, regarded as a function of the Q on the domain (5 



and is therefore continuous an d convex. As a continuous function on a compact domain, it necessarily has 
a minimum someplace on (x8), possibly on the boundary. Either e takes its minimum at a unique point or 
on some larger convex subset of (5.8). If w is a rational number one can find examples where the minimum 
is not unique. These correspond to (first-order) phase transitions where one can have two or more distinct 
ground state configurations. On the other hand, if w i s irr ational, e is a strictly convex function which 
achieves its minimum at a unique point in the simplex (5.8). The strict convexity of e in this case arises 
from the strict convexity of Gix), whose derivative with respect to x has a dense set of discontinuities on the 
interval (0, 1) when w is irrational. 

Once the I Cml are known, and thus the {zm,}, see ( 4.12| ), the ground state configuration {uj} and the hull 
function /i^, (4.3) can be obtained by means of the following construction. Define the piecewise constant 
and step periodic function 



T{z) = trn for Z„r-1 < Z < Zm, 

see Fig. ^ so that in the ground state 

Tfc = r(fcw -I- c). 

Inserting this expression in ( |3.8[ ) yields 



(5.9) 
(5.10) 



J2 9{j -k)T{ku + c)^ J2 g{l)T{]OJ + c~lu). 



(5.11) 



Comparison with (O) shows that 



E 5(0^2 -M, 



(5.12) 



^ — — oc 



where the sum yields because T has been d efined to be continuous from the left. If at each discontinuity 
one sets T{z„) equal to tm+i 



rather than t,„, (5.12) will yield 



VI. CONCLUSION 



We have shown how the problem of finding ground state configurations for a Frenkel-Kontorova model 
in a system of parabolas of identical curvature, assuming the system has a finite period, can be reduced to 
finding the minimum of a certain convex function, (5.7), over the variables {Cn} belonging to a finite simplex 
(5.8). While this second minimization problem is not trivial, the examples in our previous paper show 



8 



that the solution can be worked out in particular cases with a certain amount of effort, for instance, by 
treating the {Cn} as parameters and computing the corresponding {r]n}- An efficient numerical procedure 
for finding the minimum, which works well for N up to the order of 200, is discussed in the paper which 
follows 

It is natural to ask whether further generalizations of this type of model might prove interesting. There 
are several possibilities worth considering: First, one can imagine taking the limit N oo with P held finite 
in such a way that V converges to some smooth or piecewise smooth periodic potential, such as a cosine. 
While it is unlikely that such an approach will yield an exact solution for the ground state of an arbitrary 
smooth periodic potential, it may nonetheless provide some interesting insights, and preliminary studies of 
this problem are encouraging, see p2| . 

Second, one can imagine a limit with both N and P tending to infinity in a way which leaves their 
ratio fixed, thus yielding an example of a potential, still composed of parabolas, which is quasiperiodic or 
possesses some other deterministic form of non-periodicity. An exact solution for such a case would be quite 
interesting. 

Third, one might hope to study cases in which V consists of parabolas of two (or more) different cur- 
vatures; for example, parabolas which meet at points where their first derivatives are continuous and their 
second derivatives discontinuous. Such models have been studied numerically [^3| . An interesting analytical 
investigation of a particular case and the corresponding "tent map" has been performed by BuUett [Q| . The 
basic difficulty with non-uniform curvature is that the Green's function needed to solve the counterpart of 
(3.1) is much more difficult to calculate. Yet the fact that novel phenomena appear already when the basic 
N = 1 parabolic model is extended to = 2 suggests that any progress in this direction would be worth 
the effort. 



ACKNOWLEDGMENTS 

We are grateful to S.Aubry for bringing Ref. |2^] to our attention. One of us (R.B.G.) would like to thank 
the Institut de Hautes Etudes Scientifiques in Bures-sur-Yvette and the Service de Physique Theorique of 
the C.E.N. Saclay for their hospitality during a sabbatical year, and acknowledge financial support by the 
U. S. National Science Foundation through grant DMR-9009474. 



APPENDIX A: MINIMUM DISTANCE OF APPROACH TO DIVIDING POINTS. 

Aubry's theorems pO| for minimum energy and gro und state configurations require that V{u) be a smooth 



function, a condition violated at dividing points, (2.6), where the first derivative of the potential is discontin- 
uous. This difficulty is easily circumvented by noting that, as demonstrated below, no particle in a minimum 
energy configuration can come closer than a minimal distance > to a dividing point. Consequently one 
can always suppose, in order to apply Aubry's arguments, that V{u) has been replaced by a smooth V{u) 
identical to V{u) except for u closer to a dividing point than Sq, with V always greater than or equal to V 
(see Fig. ||). Minimum energy configurations of V are identical to those of V, as otherwise one could lower 
the energy by moving a single particle, following the reasoning given below. 

Assume that a dividing point occurs at u = 0, and that particle 1 of a minimum energy configuration falls 
in the parabola to the right of this point, at ui = 6 > 0. Ignoring an additive constant, we can write, for u 
near zero, 

with a a constant and 7 > the discontinuity in V'{u) at u = 0. 



If particles uq and U2 are held fixed, the part of the energy H, (2.1), involving ui can be written as V{ui), 
where, omitting an additive constant, 

V{u) = V{u) + u^ ~ {uo + U2)u. (A2) 
Of course V{u) must be a minimum at 
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ui = b = {uo + U2 + a)/{2 + k) > 0, (A3) 

where it takes the value 

- (uo + U2 + a))V(4 + 2K). (A4) 
If b is smaller than 7/(2 + k), V{u) also has a local minimum at 

u'l = (uo + U2 + a-7)/(2 + K) < 0, (A5) 

where it takes the value 

- (Mo + U2 + a-7)V(4 + 2K). (A6) 



Now (A6) cannot be less than (A4), as otherwise the energy would be lowered by changing ui to u[. (The 
argument is not affected by the possible prese nce of another dividing point between uq and u'l, as in that 
case the energy at u[ would be even less than (|A6|).) Consequently one has 



|uo + "2 + a - 7I = 1(2 + k)6 - 7I < uo + "2 + a = (2 + (A7) 
and thus a lower bound 

7<2(2 + k)6 (A8) 

for b. As the discontinuity 7 in V'{u) at the dividing point tj is kAIj, we conclude that no particle in a 
minimum energy configuration can come closer than 

5o = [— ^] min At, (A9) 

to one of the dividing points. 



APPENDIX B: HULL FUNCTION FOR RATIONAL WINDING NUMBER 



It is sufficient to consider the case in which P — 1 and c = 0, and as w = p/ q is fixed, the subscript lo can 
be omitted from the hull function /. As shown in |30|, a ground state configuration {«■/} is periodic, with 



for every j. Therefore, for any integer m, we can define 

fim/q) = - fc, 
where j and k are any two integers satisfying 

jp — kq ^ m. 



(Bl) 
(B2) 
(B3) 



The fact that p and q are relatively prime means that (83) always has solutions, while (Bl) ensures that 
two solutions yield the same result when inserted in (B2), and that (4^) is satisfied with P = 1: 



/(l + z) = l + /(z). 

Let s and t be two integers satisfying 

sp — tq = 1 , 

and define a second ground state configuration {vj} by means of: 

Vj = Uj+s - t. 



(B4) 
(B5) 
(B6) 
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Using @, @, and (^) one finds that: 



= / - -f mi, 

q 



/( 



TO + 1 , 

9 



TOi. 



(B7) 



A consequence of the property of "total ordering" of ground states which is proved in |Q is that if Vj < Uj 
for one value of j, this must also hold for all values of j. Applied to (B?), this means that if /[(to + l)/q] is 
less than or equ al to f{m/q) for one value of m, the same is true for all values of m, a result which clearly 
contradicts (B4). Hence 



(B8) 



for every to, and / is a strictly increasing function on the discrete set where it is defined. 
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FIG. 1. Labeling of the minima and dividing points for successive parabolas forming V{u). 



FIG. 2. The function g{x) for k = 0.5. 



FIG. 3. A hull function [z) showing the relationship of the dividing points Zn and t, 



FIG. 4. The function Q{x) for k = 0.2, w equal to the golden mean (0.618...). 



FIG. 5. The function T{z), see text. 

FIG. 6. The smooth potential t^(w), dashed curve, is identical to V(u), solid curve, except near the dividing 
points of the latter. 
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